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ABSTRACT 

In a previous paper (White and Ghosh 1998, WG98) we point out that the 
X-ray luminosity Lx of a galaxy is driven by the evolution of its X-ray binary 
population and that the profile of Lx with redshift can both serve as a diagnostic 
probe of the SFR profile and constrain evolutionary models for X-ray binaries. 

We update the WG98 work using a suite of more recently developed SFR profiles 
that span the currently plausible range. The first Chandra deep imaging results 
on L x -e volution are beginning to probe the SFR profile of bright spirals and 
the early results are consistent with predictions based on current SFR models. 
Using these new SFR profiles the resolution of the “birthrate problem” of low- 
mass X-ray binaries (LMXBs) and recycled, millisecond pulsars in terms of an 
evolving global SFR is more complete. We also discuss the possible impact of 
the variations in the SFR profile of individual galaxies. 

Subject headings: binaries: close— stars: formation— stars: evolution— galaxies: 
evolution-X-rays: galaxies-X-rays: stars 


1. INTRODUCTION 

Global star-formation rate (SFR) has undergone a strong-cosmological evolution: it was 
larger than its present value by a factor ~ 10 at z 1 , had a peak value ~ 10-100 times 
the present one in the redshift range £ ~ 1.5-3. 5, and declined again at high £ (Madau 
et al. 1996; Madau, Pozzetti k Dickinson 1998, M98; Blain, Small, Ivison k Kneib 1999, 
B99a; Blain et aZ.1999, B99b, and references therein). Details of the SFR at high redshifts 
are still somewhat uncertain, because much of the star formation at 2 z ^ 5 may be 
dust-obscured and so missed by optical surveys, but detected readily through the copious 
submillimeter emission from the dust heated by star formation (Hughes et al. 1998; Barger 
et al. 1999). 
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The X-ray emission of a normal galaxy (i.e., one without an active nucleus) is domi- 
nated by the integrated emission of the galaxy’s X-ray binary population (see e.yFabbiano 
1995). In WG98 we discussed the basic effects of an evolving SFR on the evolution of X- 
ray binary populations of galaxies, and so on that of the total X-ray emission from normal 
galaxies. The X-ray luminosities L\ of normal galaxies should show significant evolution 
(up to a factor ~ 10, depending on the LMXB evolutionary timescale), even in the relatively 
nearby redshift range z ~ 0. 5-1.0. In WG98 we also show that an evolving SFR can resolve 
the “birthrate problem” involving LMXB and their descendant “millisecond” radio pulsars 
(MRP, see Kulkarni & Narayan 1988; Lorimer 1995, L95). 

The SFR profile used in WG98— the only profile available at the time — was based on 
the optical/UV data alone. Over the past 3 years there has been considerable progress in 
our understanding of cosmic star-formation history. In addition, very deep X-ray imaging 
with Chandra is beginning to detect normal galaxies in the redshift range 2 ~ 0.5-1. 0, so 
that comparison with observations is becoming possible for the first time. In this Letter, we 
reconsider the key imprints left by the SFR evolution profile on the L_y - evolution profiles 
of normal galaxies, using the best SFR models currently available. We briefy discuss in 
this context the recent results of Brandt et al. (2001, BranOl) from the ultradeep Chandra 
imaging of the Hubble Deep Field North (HDF-N). A detailed calculation of the expected 
X-ray flux distribution of the HDF-N galaxies based on the results of this paper is reported 
by Ptak et al. (2001, P01). We also reconsider the resolution of the LMXB-MRP birthrate 
problem using the new SFR profiles. 


2. X-RAY LUMINOSITY EVOLUTION WITH EVOLVING SFR 

The total X-ray output of a normal galaxy can be modeled as the sum of those of its 
HMXB and LMXB wherein the evolution of each species “i” is described by a given timescale 
Ti (see WG98). To study the effects of the dependence of on the binary period and other 
evolutionary parameters, we run the evolutionary scheme over ranges of likely values of t, 
suggested in the literature. The (prompt) evolution of the HMXB population in response to 
an evolving star-formation rate SFR(t) is given by 

3n ™ XB(() = a h SFR(t) - " hmxbW , (1) 

dt T HMXB 

where uhmxb is the number density of HMXBs in the galaxy, and t H mxb is the HMXB 
evolution timescale. a h is a coefficient representing the rate of formation of HMXBs per 
unit SFR, given approximately by ah — | /binary /prim / sn? where /binary is the fraction of all 
stars in binaries, /£ rim is that fraction of primordial binaries which has the correct range 
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of stellar masses and orbital periods for producing HMXBs (van den Heuvel 1992, vdH92 
and the references therein), and /| N % 1 is that fraction of massive binaries which survives 
the first supernova. In the calculations reported here, we have adpoted a representative 
value thmxb ~ 5 x 10 6 yr according to current evolutionary models. Note that, in our 
introductory model here, t H mxb includes both (a) the time of evolution of the massive 
companion (~ 4 — 6 x 10 6 yr) of the neutron star from the time of the supernova (that 
produces the neutron star) to the point where the “standard” HMXB phase begins, and, (b) 
the duration ( 2.5 x 10 4 yr) of the HMXB phase (vdH92 and references therein). But since 

the second timescale is negligible compared to the first, little error is made by approximating 
this two-step process by a single-step one with an overall timescale t H mxb • 

Two basic ways of producing LMXBs have been discussed. In the cores of dense globular 
clusters, they can be produced by the tidal capture of a neutron star by a normal star. 
Over the rest of a galaxy stellar densities are insufficient for tidal capture, and LMXBs are 
produced by the evolution of primordial binaries (see, e.g., Webbink, Rappaport & Savonije 
1983; Webbink 1992). In this paper, we consider only the latter mechanism. For spiral 
galaxies, at least, this must be the dominant mechanism, since the globular-cluster LMXB 
population in such galaxies only accounts for a relatively small fraction of the total X-ray 
luminosity. 

The evolution of an LMXB from a primordial binary has two distinct stages (WG98), 
even after the supernova explosion that produces the neutron star in a post-supernova binary 
(PSNB). The PSNB first evolves on a timescale t PS nb due to nuclear evolution of the neutron 
star’s low-mass companion and/or orbital decay by gravitational radiation and magnetic 
braking, until the companion comes into Roche lobe contact and the LMXB turns on. The 
LMXB then evolves on a timescale tlmxb- Since tpsnb and tlmxb nre comparable in general, 
we have to describe the two stages separately using a formulism similar to eqn 1 (see W r G98 
for more details). 

We display evolution in terms of the redshift z, which is related to the cosmic time t by 
t 9 = 13(z + l) -3/2 , where t 9 is t in units of 10 9 yr, and a value of H 0 = 50 km s -1 Mpc -1 
has been used. We consider the suite of current SFR models detailed in Table 1 to cover a 
plausible range, using the parameterization of B99a,b. Models of the “peak” class have the 
form 

SFR peaU (^) = 2 ( 1 + exp -?-) (1 + z) p+ ^ti, (2) 

\ *max / 

while those of the “anvil” class have the form 

f (1 + z) P , Z < Zmax , 

[^ ( 1 i ** ^ /> max ■ 


SFR anv ii(^) — 


( 3 ) 
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These functional forms axe not unique, but useful, since they have a convenient low-redshift 
limit, SFR(^) oc (1 + z) p , where all SFR profiles must agree with the low -2 optical/UV 
data (M98), and the model parameters can be manipulated to mimic a wide range of star- 
formation histories (B99b). Peak-class profiles are useful for describing (a) SFRs determined 
from optical/UV observations, i.e., Madau-type profiles (M98), (called “Peak-M in Table 
1) and, (b) more general SFRs with enhanced star formation at high z , a good example of 
which is a “hierarchical” model of B99b, wherein the submillimeter emission is associated 
with galaxy mergers in an hierarchical clustering model of galaxy evolution. An\il-class 
profiles are useful for describing results of “monolithic” models. The “Gaussian” model 
(B99a,b) is an attempt at giving a good account of the SFR at both low and high 2 by 
making a composite of the Peak-G model (see Table 1) and a Gaussian starburst at a high 
redshift z p , i.e., a component 

SFR Gauss (z) = ©exp . (4) 

This is based on the IRAS luminosity function, and is devised to account for the high -2 
data, particularly the submillimeter observations (B99a). For its parameters (see Table 1), 
we have used the revised values given in B99b. In all models considered here, no galaxies 
exist for sufficiently large redshifts 2 > 10. 

Figures 1 and 2 show the evolution of HMXB and LMXB, as well as that of the total X- 
ray binary population (weighting the two components to represent the total X-ray emission 
from the galaxy). The LMXB profile peaks at redshifts ~ 1~3 later than the HMXB profile, 
which is is a charcteristic signature of SFR evolution on the X-ray binary contents of galaxies 
(WG98). At low redshifts (0 ^ 2 ^ 1), the galaxy’s X-ray emission is dominated by LMXBs, 
and at high redshifts, by HMXBs. As a result, the total Lx -profile is strongly influenced at 
high redshifts by the SFR profile, as Figure 1 shows. 

Figure 1 compares the Lx -evolution corresponding to the (Madau or peak-M) SFR 
profile and the evolutionary timescales we originally used in WG98, i.e., (a) t PS nb = 1.9 
Gyr, tlmxb =0.1 Gyr for the whole LMXB population, and'(b) t P snb = 0.9 Gyr, t L mxb = 
0.5 Gyr for the short-period systems. In Figure 2, we display the Lx-evolution for a range 
of SFR profiles— Peak-M, Hierarchical, Anvil-10, and Gaussian, the evolutionary timescales 
being held fixed at tpsnb = 1*9 Gyr, tlmxb = TO Gyr. Between them, the two figures thus 
explore the effects of (a) varying the evolutionary timescales for a fixed SFR profile, viz. 
Peak-M, and (b) varying the profile for a fixed set of evolutionary timescales. 

BranOl estimate that the average X-ray luminosity of the bright spiral galaxies at an 
average redshift 2 ss 0.5 used in their stacking analysis is about a factor of 3 higher than 
in the local Universe (2 < 0.01). This observed evolution, Lx(0.5)/Lx(0.0) ~ 3, can be 
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compaxed with the results in Table 2. The degree of evolution from z — 0 to z — 0.5 1.0 
increases from Madau-type profiles to those with additional star formation at high redshifts, 
the numbers for the Peak-M profile being in best agreement with BranOl. More sophisticated 
estimates of the expected -distribution of HDF-N galaxies, based on these evolutionary 
scenarios, axe described in P01. 


3. BIRTHRATE PROBLEM: FURTHER CONSIDERATIONS 

As discussed in WG98 the evolution of MRPs from LMXBs involves the MRP number 
density timrp, and evolution timescale tmrp ~ 3 x 10 9 -3 x 10 10 yr (Camilo et al. 1994). Our 
evolutionary scheme yields the profiles of LMXB and MRP evolution for a given SFR profile, 

giving the number ratio, N r = Wmrp /wlmxb > rate rat *°> ^ 

MRP and LMXB populations at the present epoch (z = 0). We showed in WG98 that the 
Peak-M profile yielded R. ~ 1 for the overall MRP and LMXB populations, in agreement 
with current observations (L95). However, for short-period systems (LMXB periods ^ 3 
day), this profile yielded a value R- ^ 3, smaller than the current value R- % 8 estimated 
from observations (L95). On repeating these calculations with other SFR profiles which have 
large additonal star formations at high redshifts e.g., the Gaussian and hierarchical profiles 
we now find reasonable agreement for both whole populations and short-period systems. 
These give « 6 - 8 for short-period systems, for plausible evolutionary timescales. 

A major new development in SFR research has been the first study of the star-formation 
histories of individual galaxies, and of galaxy-types. SFR profiles of individual galaxies, 
ranging from those in the Local Group to those in the HDF at redshifts 0.4 ^ z ^ 1, have 
been inferred, using several different techniques (Glazebrook etal. 1999, Abraham et al. 2000, 
Hernandez et al. 2000). For individual galaxy-types, models of spectrophotometric evolution, 
which use the synthesis code Pegase and axe constrained by deep galaxy counts (Rocca- 
Volmerange and Fioc 2000), leading to model SFR profiles. The birthrate problem was 
originally defined in terms of observations in our own galaxy and therefore we should use 
the SFR profile of our own galaxy. This profile is difficult to determine, as none of the above 
techniques is applicable to our own galaxy - but it may eventually be a tractable problem. 


4. DISCUSSION 

We have shown that different global SFR profiles within the currently admissible range 
lead to very different L_y -evolution profiles, so that the latter profiles can be an independent 
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diagnostic X-ray probe of cosmic star-formation history. The mechanism through which the 
SFR profile leaves its imprint is the interplay between the prompt evolution of HMXBs, 
which closely follows the SFR profile, and the slow evolution of LMXBs, which follow the 
SFR with a time lag of ~ 10 9 yr (WG98). Our results and a more detailed consideration by 
P01 show that the early Chandra results (BranOl) are consistent with the Peak-M profile, 
for LMXB evolution times ~ 1 Gyr. 

We can clarify the Chandra results by applying the considerations of global and in- 
dividual SFR profiles summarized in the last section. BranOl used bright spirals for their 
stacking analysis. Rocca-Volmerange and Fioc (2000) have shown that the model SFR profile 
for such (Sa-Sbc) spirals rises roughly in a Madau fashion from 2 = 0 to 2 « 1 and thereafter 
flattens to a roughly constant value ~ 12 times that at z = 0, falling again at 2 7. We 

can roughly represent this profile in the range 0 < 2 <; 7 by an anvil-type profile (see §2), 
with the parameter z ma x as given in Table 1, and the parameter p « 2.7. On running our 
evolutionary scheme for such a profile with the timescales t PS nb = 1-9 Gyr, tlmxb = 1-0 
Gyr, as in Figure 2, we get L,y(0.5)/L.y(0.0) = 3.3, and L,y(10)/Z/_y(0.0) = 5.4. This is 
in good agreement with both the BranOl results and the Peak-M results summarized in §2, 
and so provides a natural explanation for why the Peak-M profile appears to give a good 
account of the BranOl results. In effect, the BranOl analysis may be probing the SFR profile 
of bright spirals in HDF-N. While the global analysis presented here is conceptually the most 
appropriate statistical descrption of a field like the HDF, we must be able to sample all types 
of galaxies in order to probe the global SFR observationally. 

The X-ray probe considered here is a “fossil record” of SFR, as LMXBs preserve this 
record on a timescales ~ 1 Gyr. Conversely, as proving grounds for theories of X-ray binary 
evolution, this probe enjoyes a unique position, since such theories have been built almost 
wholly on the basis of our experience in the current epoch. This new probe provides a new 
means of watching LMXB evolution unfold over cosmic time. 

It is a pleasure to thank L. Angelini, R. Griffiths, R. Mushotzky, and A. Ptak for 
stimulating discussions. 
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Table 1: Stax Formation Rate (SFR) Profiles® 


Model 

& max 

V 

Comments 

Peak-M 

0.39 

4.6 

Madau profile 

Hierarchical 

0.73 

4.8 

Hierarchical clustering model 6 

Anvil-10 

1.49 

3.8 

Monolithic models 

Peak-G 

0.63 

3.9 

“Peak” part of composite Gaussian 

Gaussian 

N/A 

N/A 

Gaussian starburst® added at high-x 


“Model parameters taken from B99a,b. 

b We have chosen the B99b model with a dust temperature 45 K. Note that the parameter p 
for the SFR equals 3/2 plus the value of the parameter p occurring in equation (16) for the 
merger efficiency in B99b. 

c Parameters of Gaussian starburst component (see eq.[6] of text), are from the modified 
model given in B99b, i.e., z p = 1.7, a — 1.0 Gyr, and © = 70. 


Table 2: Evolution of X-ray Luminosity L x 


Model 

tpsnb 

tlmxb 

Lxi 0-5) 

Xx(O.O) 

Lx( 10) 
Lx (0.0) 

Peak-M 

1.9 

0.1 

3.9 

5.4 

Peak-M 

0.9 

0.5 

4.6 

6.8 

Peak-M 

1.9 

1.0 

3.4 

4.1 

Hierarchical 

1.9 

1.0 

6.2 

11.3 

Anvil-10 

1.9 

1.0 

5.4 

10.1 

Gaussian 

1.9 

1.0 

7.5 

16.0 
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Fig. 1. — Evolution of HMXB population (dotted line), LMXB population (dash-dotted line), 
and the total X-ray luminosity Lx (solid line) of a galaxy with a given SFR (dashed line). As 
absolute ordinate scales are irrelevant for these evolutionary profiles, they have been adjusted 
for convenience of display: Lx always starts at 0.001 at z = 0, so that its evolution can be 
immediately read off the figure, and SFR always starts at 0.01 at z = 0, so that different SFR 
profiles can be readily compared. This figure is for the Peak-M profile (see text), showing 
the effect of varying the evolutionary timescales tpsnb and tlmxb, whose values are written 
on each panel. The timescales used here are those used in WG98, corresponding to whole 
LMXB populations and short-period systems. 

Fig. 2. — S am e as Figure 1, but showing the effect of varying the SFR profile. The evolu- 
tionary timescales are kept fixed at tpsnb = 1-9 Gyr and tlmxb = 1-0 Gyr for all cases, and 
SFR profiles from Table 1 are used. Each panel is labeled by the name of its SFR profile. 
Evolutionary factors from this and the previous figure are collected in Table 2 and described 
in the text. 
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